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The question of a modification of the running gauge coupling of (non-) abelian gauge theories by 
an incorporation of the quantum gravity contribution has recently attracted considerable interest. 
In this letter we perform an involved diagrammatical calculation in the full Einstein- Yang-Mills 
system both in cut-off and dimensional regularization at one loop order. It is found that all gravita- 
tional quadratic divergencies cancel in cut-off regularization and are trivially absent in dimensional 
regularization so that there is no alteration to asymptotic freedom at high energies. The logarithmic 
divergencies give rise to an extended effective Einstein- Yang-Mills Lagrangian with a counterterm 
of dimension six. In the pure Yang-Mills sector this counterterm can be removed by a nonlin- 
ear field redefinition of the gauge potential, reproducing a classical result of Deser, Tsao and van 
Nieuwenhuizen obtained in the background field method with dimensional regularization. 
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Pcrturbativcly quantized general relativity is well 
known to be a non-renormalizablc theory due to the neg- 
ative mass dimension of its coupling constant k [i[ . The 
coupling to matter fields does generically not improve 
the situation [3, 0, Q j although the possibility of pertur- 
bative finiteness of the maximally supersymmetric grav- 
ity theory is still open Hence, Einstein's theory of 
gravity does not constitute a fundamental theory as its 
non-renormalizability necessitates the inclusion of an in- 
finite set of higher dimension counterterms in the per- 
turbative quantization process. Nevertheless, as advo- 
cated by Donoghue it may be treated as an approx- 
imation to a fundamental theory of quantum gravity by 
the methodology of effective field theories in order to de- 
scribe interactions at scales well below the Planck mass 
Mp ~ 1/k ~ IQi^GeV. 

In an interesting recent paper Robinson and Wilczek 
reported on a calculation in Einstein gravity coupled 
to Yang-Mills theory, where the quantum gravity con- 
tributions to the running of the Yang-Mills coupling g 
were investigated at one-loop order. Interestingly a neg- 
ative gravitational contribution to the Callan-Symanzik 
/3 function, which quantifies the flow of the Yang-Mills 
coupling with the energy scale E, was found 
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with oo = —3/2 in our conventions for the gravitational 
coupling K^, see ([2]). Irrespective of the non-gravitational 
value of bo this would render any gauge theory assymp- 
totically free at energies E close to the Planck mass (in- 
cluding e.g. pure U{1) Maxwell theory). Responsible for 
this effect were quadratic divergencies in one-loop graphs 
containing the graviton propagator 0. As a number of 
discussed scenarios of physics beyond the standard model 
contain higher dimensional gravity theories with small 
gravitational scales, such an effect might be observable 
at the Large Hadron Collider, as was qualitatively stud- 
ied in M- 



However, two recent works have cast doubts on the 
results of by reconsidering this effect in Einstein- 
Maxwell theory. The authors of @ used the back- 
ground field method, which contains subtle dependen- 
cies on gauge conditions leading to gauge dependent re- 
sults of the effective action. To our undestanding the 
essence of the problem in Q lies in expanding the grav- 
itational field about flat Minkowski space and the gauge 
field about a nontrivial solution of the Hat space Yang- 
Mills field equations, which docs not constitute a solution 
of the coupled Einstein- Yang-Mills system. Indeed in Q 
Pietrykowski repeated the analysis of @ in the abelian 
theory employing a parameter dependent gauge condition 
and obtained a gauge dependent result. In a very recent 
analysis Toms [10[ pointed out the aforementioned sub- 
tleties in the background field approach and studied the 
problem in Einstein-Maxwell theory employing an elab- 
orated gauge invariant and gauge condition independent 
formulation of the background field method due to Vilko- 
visky and De Witt [ll[ . Here dimensional regularization 
was used and a vanishing gravitational contribution to 
the gauge coupling /3-function was found. 

In view of these discussions it should be stressed that 
the question of one loop divergencies in the Einstein- 
Yang-Mills theory was already studied more than 30 
years ago in the classical work of Deser, Tsao and van 
Nieuwenhuizen In a background field approach con- 
sidering fluctuations about a general gravitational and 
gauge field background the necessary one loop countert- 
erms were computed by these authors in dimensional 
regularization. In the pure Yang-Mills sector only one 
dimension-six counterterm of the form Tr[(Z)^F^'^)^] 
arose, which moreover must be removed through a non- 
linear field redefinition of the gauge field. Hence already 
these results were in contradiction with the findings of 
0. However, one might be worried about the use of 
dimensional regularization in all these computations, it 
being insensitive to quadratic divergencies which at the 
same time give rise to a nonvanishing ag in the Callan- 
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Symanzik /3-function ([T]) as reported in Q using a cutoff 
prescription. 

In view of the recent interest in the problem, its poten- 
tial experimental relevance and the discussed technical 
controveries we have conducted a conceptually straight- 
forward but involved diagrammatical calculation in the 
full Einstein- Yang-Mills system in cut-off and dimen- 
sional regularization. For this the gravitational contribu- 
tions to the gluon self-energy and vertex function were 
evaluated at one loop order. The relevant diagrams are 
listed in figure [TJ We find that all quadratic divergen- 
cies cancel in cut-off regularization, and they are triv- 
ially absent in dimensional regularization. The logarith- 
mic divergent contributions are substracted by precisely 
the dimension-six countcrtem found by Descr, Tsao and 
van Nicuwcnhuizen Q confirming their results and show- 
ing explicitly the non-renormalizability of the Einstein- 
Yang-Mills theory. In particular the /3-function of the 
Yang-Mills coupling constant g receives no gravitational 
contributions at the one loop order, in contradiction to 
the findings of 0], i.e. oq = in H]). 

The starting point of our analysis is the Einstein- Yang- 
Mills theory 
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and the two gluon-two graviton vertex 
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where we have defined pf^''-"'^ = l[n^'°''q^P + r|^'P■q'''^ 
jy^'^yy"''). Using these expressions we have computed the 
divergent pieces of the two-gluon graphs (a) and (b) of 
figure [T] with an incoming momentum of q to be 
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with the Ricci scalar R and the Yang-Mills field strength 
F^u = dfj^Ai, — dvA^ ~i 9 ^ -^A ■ In the renormalization 
process we will be led to add additional dimension-six 
operators to this theory, as discussed above. The metric 
tensor is split into a flat Minkowski background jy^^ = 
diag(l, —1, —1, —1) and the graviton field /i^j/ 



(3) 



The Einstein- Yang-Mills lagrangian is then expanded up 
to second order in h^i, in order to read off the relevant 
propagator (modulo gauge fixing) and gluon-graviton 
vertices occuring in figure [T] We work in Feynman gauge 
for the gluons and in harmonic (de Donder) gauge for the 
gravitons with the gauge condition d^h^j^^, 
resulting graviton propagator reads 
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where we quote both the cut off {\k'^\ < A^) and di- 
mensional (d — A — e) regularized results. For details 
of the calculation see Note the cancellation of the 

quadratic divergence which already leads to the absence 
of gravitational corrections to the /? function of g in the 
abelian theory. 

For an abelian gauge theory this would be all there is 
to do, as there are no cubic gauge field vertices. In order 
to substract the remaining divergence in ([7]) at a given 
renormalization point /i we have to augment the Einstein- 
Yang-Mills lagrangian ^ by novel dimension-six coun- 
tertems. Taking into account the Bianchi identity there 
are naively three possible structures 



Oi = Tr[ (D^F^pf ] , O2 = Tr[ {D^F'^^f ] , 
03 = i TriF/F.^F/]. (8) 



However, it turns out that they are linearly related up to 
total derivative terms 



with 7^'''"'' = ^{ififf^ + Tj'^'^if"'). As gravitational 
ghosts do not couple to gluons they do not appear in 
our one-loop computations. We furthermore note the 
two gluon-one graviton vertex 
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■ Oi ~ 2 g O3 + total derivatives . 
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We are thus led to add the terms di Oi and d2 O2 to 
our lagrangian, where di^2 are the corresponding cou- 
pling constants. Note that the term d2 O2 in the ex- 
tended effective lagrangian is proportional to the low- 
est order equations of motion D^j^F^^^ = and could 
be removed by a field redcfiniton of the gauge field 
— > — d2 D^F'^ ^/2 up to higher dimension oper- 
ators. 
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FIG. 1: Graviton loop corrections to the gluon two and three point functions. 



The new two gluon vertices of di Oi and d2 O2 are 



tributions in cut-off and dimensional regularization 
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The associated counterterms to these vertices can be used 
to cancel the logarithmic poles of (a) + (&) in ([7]). How- 
ever, due to the identical tensor structure of pU|) and 
(jlip these two counterterms are only fixed up to one free 
parameter. 

In order to fix this free parameter we move on to the 
thrcc-gluon graphs (c) , (d) and (c) which probe the non- 
abelian sector of the theory. Expanding Ceym to cubic 
order in gluons and up to quadratic order in gravitons 
we read off the relevant three gluon-one graviton vertex 
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+T]^"''qP'^'^{p- qf^ 
+cycl(Ai,p; V, q; p,k) 



and the involved three gluon-two graviton vertex 



+ {(«,/?)<- (7, <S)} 
+cycl(pt,_p; i^,q;p,k)] . 

(13) 

With these the evaluation of the three-gluon graphs of 
figure 1 can be performed, and wc find the divergent con- 
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where the external gluons carry the labels {fJ.,p', v, q; p, k) 
and the dots refer to a symmetrization in these labelings. 
Summing these contributions up we again observe a can- 
cellation of the quadratic divergencies arising in cut-off 
regularization. In order to subtract the remaining log- 
arithmic divergencies in the three-gluon amplitudes we 
need to consider the three point vertices of di Oi and 
d2 O2 emerging from ([5]) 

vh 

^digr'"-i,npP{4p-q + 2p-k) 
p " ~qP{4q-p + 2q-k)) ' 

k -2 (fc^r(p-g)'' + ...)!, (17) 



= d2gr'^[v^''{pP{2p -q + p-k + Zq-k) 
-qP{2q-p + q- k + 3p- k)) + ... 
(02)^ /-ia -{k^'k''{p-q)P + ...) 
P -SipPq^k" -p^qPkP)] . 

^ (18) 
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Remarkably it turns out that O2 alone provides the right 
tensor structures to remove all the divergences of p^ - 
(fT6|l . Renormalization of the two and three point gluon 
functions is then performed by considering the extended 
effective Einstein- Yang-Mills theory augmented by 3,2 O2 
plus all the associated counter terms. The complete ef- 
fective lagrangian at the one-loop level is then of the 
schematic form 



C =£oxt + CcT , 

-^CeYM, rcn + ^2 Tr[(i?^^^^'')2] , 

CT =62 {dAf + g <5?» A'dA + ^i^^ {d" Af 
+ gSl'^^ d'^AdAA + 0{A^). 
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The computed divergencies are then substracted at a 
renormalization scale /i through a choice of renormaliza- 
tion conditions. This leads to the following dependen- 
cies of the dimension-six countertcrms 
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where the first equality S^"^^ \o{k.^) = '^i'^^ \o{k^) arises as a 
necessary condition on the renormalization of the gluon 
selfenergy and vertex function and constitutes an inde- 
pendent consistency check of our results as required by 
the Slavnov- Taylor- Ward identity. Also note that this 
countertem is in prescise agreement to the one found by 
Deser, Tsao and van Nieuwenhuizen Q with the back- 
ground field method. In particular the Yang-Mills vertex 
counterterm Si and wavefunction renormalization S2 re- 
ceive no contributions at order k^. The countertcrms 
(PO)) then lead to the following relations of renormalized 
g, d2 to bare g^^ d2.o couplings at one loop order 
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This then provides the following gravitational contribu- 
tions to the P functions of the gauge coupling g and the 
novel coupling d2 
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However, the renormalization of d2 through gravitational 
interactions is not of particular relevance as this coupling 
can be removed through a nonlinear field redefinition as 
mentioned above. 

Let us qualitatively extend our discussion to the gluon 
four point function. A simple dimensional considera- 
tion of loop integrals for the gluon four point function 
yields a new operator structure K^Tr(i^^), which can- 
not be removed by the above quoted field redefinition. 
This is in agreement with the findings of where it 
was shown that such terms appear as the square of the 
energy-momentum tensor, (r^^)^. 

Interesting further extension of this work would be to 
consider large extra dimension scenarios, where a cut- 
off prescription appears mandatory, as well as the cou- 
pling to additional matter fields. The case of gauged 
N extended supergravities was considered in [l3| in the 
background field method with dimensional regulariza- 
tion, where a vanishing of the beta function for A'^ > 4 
was found. 
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